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Squared eigenfunctions are quadratic combinations of Jost functions and adjoint Jost functions 
which satisfy the linearized equation of an integrable equation. They are needed for various studies 
related to integrable equations, such as the development of its soliton perturbation theory. In this 
article, squared eigenfunctions are derived for the Sasa-Satsuma equation whose spectral operator 
is a 3 x 3 system, while its linearized operator is a 2 x 2 system. It is shown that these squared 
' eigenfunctions are sums of two terms, where each term is a product of a Jost function and an 

^ , adjoint Jost function. The procedure of this derivation consists of two steps: first is to calculate the 

■ variations of the potentials via variations of the scattering data by the Riemann-Hilbert method. The 

' second one is to calculate the variations of the scattering data via the variations of the potentials 

through elementary calculations. While this procedure has been used before on other integrable 
equations, it is shown here, for the first time, that for a general integrable equation, the functions 
appearing in these variation relations are precisely the squared eigenfunctions and adjoint squared 
eigenfunctions satisfying respectively the linearized equation and the adjoint linearized equation of 
£ — . the integrable system. This proof clarifies this procedure and provides a unified explanation for 

previous results of squared eigenfunctions on individual integrable equations. This procedure uses 
primarily the spectral operator of the Lax pair. Thus two equations in the same integrable hierarchy 
, will share the same squared eigenfunctions (except for a time-dependent factor). In the Appendix, 

the squared eigenfunctions are presented for the Manakov equations whose spectral operator is 
closely related to that of the Sasa-Satsuma equation. 



> 
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I. INTRODUCTION 



Squared eigenfunctions are quadratic combinations of Jost functions and adjoint Jost functions which satisfy the 
linearized equation of an integrable system. This name was derived from the fact that, for many familiar integrable 
. equations such as the KdV and NLS equations, solutions of the linearized equations are squares of Jost functions of 
7— I ' the Lax pairs [3, Sill, 0,0. Squared eigenfunctions arc intimately related to the integrable-equation theory. For 
\ instance, they are eigenfunctions of the recursio n op erator of integrable equations [si. I. llCll. E|. They also appear 
as self-consistent sources of integrable equations jl3l . Ej . A more important application of squared eigenfunctions is 
^\ \ in the direct soliton perturbation theory, where squared eigenfunctions and their closure relation play a fundamental 
O ■ role 

Squared eigenfunctions have been derived for a number of integrable equations including the KdV hierarchy, the 
• AKNS hierarchy, the derivative NLS hierarchy, the sine-Gordon equation, the massive-Thirring model, the Benjamin- 

' Ono equation, the matrix NLS equations, the Kadomtsev-Petviashvili equation, etc. 1 111, [H, El EH, El EH- 
Several techniques have been used in those derivations. One is to notice the commutability relation between the 
linearized operator and the recursion operator of the integrable equation [ToL El fl9| . Thus squared eigenfunctions 
are simply eigenfunctions of the recursion operator. The drawback of this method is that one has to first derive the 
recursion operator of the integrable equation and determine its eigenfunctions, which can be highly non-trivial for 
many equations (see [2(| for instance). Another technique is to first calculate variations of the scattering data via 
the variations of the potentials, then use the Wronskian relations to invert this relation to get the variations of the 
potentials via variations of the scattering data 0]. This Wronskian-relation technique involves some ingenious steps 
and could run into difficulties too if the problem is sufficiently complicated 0]. The third technique is related to the 
second one, except that one directly calculates the variations of the potentials via variations of the scattering data 
by the Riemann-Hilbert method ll]. This third method is general and conceptually simpler. Regarding the second 
and third methods however, one question which was never clarified is why the functions appearing in the variation 
relations are indeed squared eigenfunctions and adjoint squared eigenfunctions which satisfy respectively the linearized 
equation and adjoint linearized equation of an integrable system. In all known examples, this was always found to be 
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true by direct verifications. But whether and why it would remain true for the general case was not known. 

Recently, we were interested in the Sasa-Satsuma equation which is relevant for the propagation of ultra-short 
optical pulses [U, [22|, [23[ . This equation is integrable [23| ■ Two interesting features about this equation arc that 
its solitons are embedded inside the continuous spectrum of the equation (la . |24| , and their shapes can be double- 
humped for a wide range of soliton parameters [23 1. One wonders how these double- humped solitons evolve when the 
equation is perturbed — a question which is interesting and significant from both physical and mathematical points 
of view. This question can be studied by a soliton perturbation theory. Due to the embedded nature of these solitons, 
external perturbations will generally excite continuous-wave radiation which is in resonance with the soliton. So a key 
component in the soliton perturbation theory would be to calculate this continuous- wave radiation, for which squared 
eigenfunctions are needed (a similar situation occurs in the soliton perturbation theory for the Hirota equation [H|). 
A special feature of the Sasa-Satsuma equation is that, while the linearized operator of this equation is 2 x 2, its 
spectral operator is 3 x 3. How to build two-component squared eigenfunctions from three-component Jost functions 
is an interesting and curious question. As we shall see, only certain components of the 3x3 Jost functions and their 
adjoints are used to construct the squared eigenfunctions and the adjoint squared eigenfunctions of the Sasa-Satsuma 
equation. To calculate these squared eigenfunctions for the Sasa-Satsuma equation, we have tried to directly use 
the first method mentioned above. The recursion operator for the Sasa-Satsuma equation has been derived recently 
[20| . However, that operator is quite complicated, thus its eigenfunctions are difficult to obtain. The second method 
mentioned above is met with difficulties too. Thus we choose the third method for the Sasa-Satsuma equation and 
demonstrate here certain advantages of this method. 

In this paper, we further develop the third method, using the Sasa-Satsuma equation as an example. We show that 
the functions appearing in the expansion of the variations of the potentials arc always the squared eigenfunctions 
which satisfy the linearized equation of an integrable system, and that the functions appearing in the formulae for 
the variations of scattering data are always the adjoint squared eigenfunctions which satisfy the adjoint linearized 
equation of an integrable system. In addition, given these two relations between the variations of the potentials and 
the variations of scattering data, there naturally follows the closure relation for the squared eigenfunctions and their 
adjoints, as well as all the inner-product relations between the squared eigenfunctions and their adjoints. Thus no 
longer is it necessary to grind away at calculating these inner products from the asymptotics of the Jost functions. 
Rather one can just read off the values of the nonzero inner products from these two variation relations. This clarifies 
the long-standing question regarding squared eigenfunctions in connection with the linearized integrable equation, 
and streamlines the third method as a general and conceptually-simple procedure for the derivation of squared 
eigenfunctions. We apply this method to the Sasa-Satsuma equation, and find that it readily gives the squared 
eigenfunctions and adjoint squared eigenfunctions. The squared eigenfunctions for the Sasa-Satsuma equation are 
sums of two terms, with each term being a product of a component of a Jost function and a component of an 
adjoint Jost function. This two-term structure of squared eigenfunctions is caused by the symmetry properties of 
scattering data of the Sasa-Satsuma equation. It should be noted that our derivation uses almost exclusively the 
spectral operator of the Lax pair, thus all integrable equations with the same spectral operator (such as the Sasa- 
Satsuma hierarchy) will share the same set of squared eigenfunctions as we derived here (except for a time-dependent 
factor which is equation-specific). An additional benefit of the method we used is that for two integrable equations 
with similar spectral operators, the derivation of their squared eigenfunctions will be essentially the same. Thus one 
can get squared eigenfunctions for one equation by minor modifications for the other equation. As an example, we 
demonstrate in the Appendix how squared eigenfunctions for the Manakov equations can be easily obtained by minor 
modifications of our calculations for the Sasa-Satsuma equation. 



II. THE RIEMANN-HILBERT PROBLEM 



To start our analysis, we first formulate the Ricmann-Hilbcrt problem for the Sasa-Satsuma equation which will be 
needed for later calculations. The Sasa-Satsuma equation is 

ut + Uxxx + 6\u\ 2 u x + 3u(\u\ 2 ) x = 0. (1) 

Its spectral (scattering) problem of the Lax pair is [23[ 

Y x = -iQKY + QY, (2) 

where Y is a matrix function, A = diag(l, 1, —1), 

/ u \ 

Q = u* ) (3) 
\ -u* -u / 
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is the potential matrix, the superscript "*" represents complex conjugation, and ( is a spectral parameter. In this 
paper, we always assume that the potential u{x) decays to zero sufficiently fast as x — ► ±00. Notice that this matrix 
Q has two symmetry properties. One is that it is anti-Hcrmitian, i.e. Qt = —Q, where the superscript "|" represents 
the Hermitian of a matrix. The other one is that 

/0 1 0\ 

aQa = Q*, where a = cr" 1 = 1 . (4) 

V° 1/ 

Introducing new variables 

J = YE-\ E = e- lCAx , (5) 

then Eq. (J2J) becomes 

J* = -<C[A, J] + QJ, (6) 

where [A, J] = A J — J A. The matrix Jost solutions J±(x, £) of Eq. ([6]) are defined by the asymptotics 

J± — ► 1, as 2; — > ±00, (7) 

where "1" is the unit matrix. Here the subscripts in J± refer to which end of the x-axis the boundary conditions are 
set. Since trQ = 0, using Abel's formula on Eq. ((6|), we see that det J± = 1 for all x. In addition, since J±E are 
both solutions of the linear equations @ , they are not independent and are linearly related by the scattering matrix 
S(C) = MC)]: 

J- = J+ESE- 1 , dctS = 1 (8) 

for real £, i.e. ( 6 1 (for non-real C, certain elements in S may not be well defined). 

Due to the two symmetry properties of the potential matrix Q, the Jost solutions J± and the scattering matrix S 
satisfy the corresponding symmetry properties. One is the involutive property. Since the potential Q is anti-Hcrmitian, 
we see that J±(C*) and ^± X (C) both satisfy the adjoint equation of j6]). In addition, we see from Eq. ([7]) that J±((*) 
and J± 1 (C) have the same large- a; asymptotics, thus they are equal to each other: 

4(D = j ± 1 (c). (9) 

Then in view of Eq. ([5]) , we see that 

s\e) = s- 1 ((). (10) 

To derive the other symmetry properties of the Jost solutions and the scattering matrix, we notice that due to the 
symmetry (|4]) of the potential Q, it is easy to see that aJ±(— (*)<r also satisfies Eq. (|6|). Then in view of the 
asymptotics (j7]l . we find that Jost solutions possesses the following additional symmetry 

J±(C)=a4(-C>. (11) 
From this symmetry and the relation l|8|). we see that the scattering matrix S possesses the additional symmetry 

s(o = *s*(-cy (12) 

Analytical properties of the Jost solutions play a fundamental role in the Riemann-Hilbert formulation for the 
scattering problem ([6]). For convenience, we express J± as: 

J_ = $£-\ (13) 
j+ = ^e-\ * = hk,^,Va], (14) 

where $ and $ are fundamental solutions of the original spectral problem ([2]) and are related by the scattering matrix 
as 



(15) 
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Then in view of the boundary conditions l[7|). the spectral equation ([6]) for J± can be rewritten into the following 
Volterra-type integral equations: 

J_(C;z) =t + T e* K{ y-^Q{y)J^ y y^-yUy, (16) 



/■ oc 

J + (C; x)=l- e^-^Qiy) J_(C; y)e< A ^-^dy. 



(17) 



These integral equations always have solutions when the integrals on their right hand sides converge. Due to the 
structure ((3]) of the potential Q, we easily see that Eq. (fT6|) for the first and second columns of J- contain only the 
exponential factor e l ^ x ~ v ^ which decays when ( is in the upper half plane C+, and Eq. (fT71) for the third column 
of J + contains only the exponential factor e l ^ v ~ x ^ which also falls off for £ £ C+. Thus these three columns can be 
analytically extended to ( £ C+. In other words, Jost solutions 



are analytic in £ £ C+, where 



P+ = [</)!, 02, ip 3 ]e i€Ax = J-Hi + J+H 2 



Hj. = diag(l,l,0), F 2 = diag(0,0,l) 



(18) 



(19) 



Here the subscript in P refers to which half plane the functions are analytic in. From the Volterra integral equations 
for P + , we see that 



P + (x,Q->t as CeC+^oo. 
Similarly, Jost functions [?/>i, ip2, <j)3]e l ^ Ax are analytic in C £ C_, and their large-C asymptotics is 

[?Puip2,fo]e lCAx -> 1 as (eC-^oo. 
To obtain the analytic counterpart of P+ in C_, we consider the adjoint spectral equation of ([£ 

K x = -i([A,K] - KQ. 
The inverse matrices J^ 1 satisfy this adjoint equation. Notice that 

JZ 1 = P*" 1 , JZ 1 = E^- 1 . 
Let us express and collection of rows, 



EE * = 



1p3 



(20) 
(21) 
(22) 
(23) 

(24) 



where the overbar refers to the adjoint quantity, then by similar techniques as used above, we can show that the first 
and second rows of JZ 1 and the third row of JZ 1 arc analytic in ( £ C_, i.e. adjoint Jost solutions 



P_ 



02 
1p3 



HlJZ 1 + H 2 JZ 1 



are analytic in ( £ C_. In addition, their large-C asymptotics is 

P_(a;,C)-»l as C e C_ -> oo. 



(25) 



(26) 



Similarly, the first and second rows of J, 1 and the third row of J_ , i.e. e ^"ipi, e ^V^j and e^fo, are analytic 
in £ £ C+, and their large-£ asymptotics is 



^2 
4>3 



1 as C £ 



(27) 
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In view of the involution properties ((9|) of J± , we see that the analytic solutions P± satisfy the involutive property as 
well: 



pUC) = P-((). 



(28) 



This property can be taken as a definition of the analytic function P_ from the known analytic function P + . 

The analytic properties of Jost functions described above have immediate implications on the analytic properties 
of the scattering matrix S. Let us denote 



S- 1 (C) = 5(C) = [5«(C)], 



then since 



S = = 



02 
03 



■>1,<P2,<P3 



S = = 



[^1,^2,^3], 



(29) 



(30) 



we see immediately that s u , s 12 , s 21 , s 22 and s 33 can be analytically extended to the upper half plane ( G C+, while 
sn, S12, S21, S22 and S33 can be analytically extended to the lower half plane £ G C_. In addition, their large-£ 
asymptotics are 



and 



Sll S12 
S21 s 22 



Sll S12 

S21 s 2 2 



1, s 33 ->• 1, as C e 



1, s 33 1, as C e 



00, 



(31) 



(32) 



Hence we have constructed two matrix functions P + and P_ which are analytic in C + and C_ respectively. On the 
real line, using Eqs. ©, (fTS|) . I|25p. and l]29p. we easily see that 



where 



P-(C)P+(C) = G(0, CeR, 



1 S13 

G = P(Pi +H 2 S){H l + S- 1 H 2 )E- 1 = E | 1 s 23 | E' 1 . 

s 31 s 32 1 



(33) 



(34) 



Equation (|33[) determines a matrix Ricmann-Hilbcrt problem. The normalization condition for this Ricmann-Hilbcrt 
problem can be seen from (|2"0|) and (f2T))) as 



■PiC^O-^l as (e 



(35) 



If this problem can be solved, then the potential Q can be reconstructed from an asymptotic expansion of its solution 
for large £. Indeed, writing P + as 



p+( x ,o = 1 + c _1 -p (1) (a;) + r 2 ^p (2) (^) + ^(r 3 ), 

inserting it into Eq. ([6]) and comparing terms of the same order in we find at O(l) that 

Q = *[A,P«]. 

Thus the potential Q can be reconstructed from pW. At 0(C _1 ), we find that 

P« =-i[AP«] + Qp( 1 ). 



(36) 
(37) 
(38) 



From the above two equations as well as the large- x asymptotics of P+(x, £) from Eq. ([7]), we see that the full matrix 
P« is 



K?/)| 2 */ J^ oo u 2 (y)dy u(x) 
Il oc u* 2 (y)dy J^luiy^dy u*(x) 

u*(x) u(x) 2f™\u(y)\ 2 dy 



(39) 
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This matrix, together with Eq. (|36|) . gives the leading-order asymptotic expansions for the analytical functions P+. 
Leading-order asymptotic expansions for P_ can be similarly derived. 

Elements of the scattering matrix S arc not all independent. Indeed, since the potential matrix Q contains only one 
independent function u(x), matrix S should contain only one independent element as well. To get the dependence of 
scattering coefficients, we first notice from equation SS = 1 that 

«33(C)a33(C) = l-*3i(0*i3(0-*3a(C)*23(C) > C G R. (40) 

Since S33 and S33 are analytic in C_ and C + respectively, the above equation defines a (scalar) Ricmann-Hilbcrt 
problem for (S33, S33), under the canonical normalization conditions which can be seen from (|31[) - (|32[) . Thus (S33, S33) 
can be uniquely determined by the locations of their zeros in C± together with the scattering data (S31, S32, S13, S23) 
on the real line £ G R by the Plcmelj formula. Similarly, from equation SS = 1 we have 

SU S 12 \ (SU. S12 ^ _ t _ / SIS \ (S31) S32); c g Rj (41) 



S21 S22 j \ S21 s 22 / V s23 . 

which defines another (matrix) Riemann-Hilbcrt problem for (sn, S12, S21, S22) and (sn, S12, S21, s 22 ), under the 
canonical normalization conditions which can be seen from ([3~T|) - ([32|) . Thus these analytical scattering elements 
can be determined by the scattering data (S31, S32, S13, S23) on the real line ( £ R as well. The scattering data 
( s 3ii S32, S13, S23) themselves are dependent on each other by the symmetry relations (|I0[) and (112p . Specifically, for 
£ £ R, 513(C) = s 3i(£); 523(C) = s 32(£)i an d S32(0 = s 3i(— 0- Thus the scattering matrix S indeed contains only a 
single independent element as we would expect. 

In the next section, we will derive squared eigenfunctions for the Sasa-Satsuma equation. For the convenience of 
that derivation, we introduce the following notations: 

S31 S32 _ S13 _ S 2 3 ,, n , 

pi = — , P2 = — , pi = —, p 2 = —. (42) 

s 33 s 33 s 33 s 33 

Due to the symmetry conditions (|10[) and (|12p of the scattering matrix 5", we see that pk and p^ satisfy the symmetry 
conditions 

p fc (C) = Pfc(C),fc = i,2, pi(C) = P2(-C), p 2 (C)=p 1 (-C), CeR. (43) 

In the next section, symmetry properties of $ and ^ will also be needed. Due to the symmetry conditions ([9]) and 
(fTTj) of the Jost solutions, $ and ^> satisfy the symmetry conditions 

$t(o = * _1 (o. * t (a = *- 1 (o. (44) 
* ± (0 = ff*i(-c> ) *±(o = ff*±(-c*)o-. (45) 



and 



Symmetry (]44[) means that 

4(0 = ^(0, 4(0=^(0, fc = l,2,3, (46) 
while symmetry (j45]l. together with (|46|) . means that 

^i(0 = ^(-C), ^2(C) = ^(-C), ^3(0 = ^3 (-0- (47) 

Here the superscript "T" represents the transpose of a matrix. These symmetry properties will be important for 
deriving the final expressions of squared eigenfunctions for the Sasa-Satsuma equation in the next section. 



III. SQUARED EIGENFUNCTIONS AND THEIR CLOSURE RELATION 

In this section, we calculate the variation of the potential via variations of the scattering data, then calculate 
variations of the scattering data via the variation of the potential. The first step will yield squared eigenfunctions, 
and it will be done by the Riemann-Hilbert method. The second step will yield adjoint squared eigenfunctions, and 
it will be done using basic relations of the spectral problem ([2]). For the ease of presentation, we first assume that S33 
and S33 have no zeros in their respective planes of analyticity, i.e. the spectral problem (|2|) has no discrete eigenvalues. 
This facilitates the derivation of squared eigenfunctions. The results for the general case of S33 and S33 having zeros 
will be given in the end of this section. 
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A. Variation of the potential and squared eigenfunctions 



In this subsection, we derive the variation of the potential via variations of the scattering data, which will readily 
yield the squared eigenfunctions satisfying the linearized Sasa-Satsuma equation. Our derivation will be based on the 
Ricmann-Hilbcrt method. 

To proceed, wc define the following matrix functions 



(48) 



F + = P+ diag(l, 1, F_ = PI 1 diag(l, 1, s 33 ). 

The reason to introduce diagonal matrices with S33 and S33 in F± is to obtain a new Ricmann-Hilbcrt problem (|49[) 
with a connection matrix G which depends on pj~ and pk rather than Sjj and Sij. This way, the variation of the 
potential will be expressed in terms of variations in pk and p^. When S33 and S33 have no zeros in C+ and C_ 
respectively, then F± as well as FjT 1 are analytic in C±. On the real line, they are related by 



F + (C) = F_(OG(C), (e 



where 



1 



1 



1 pi 



G = diag(l,l, — )(?diag(l,l,— )=E 1 p 2 



E~ 



S33 



S33 



(49) 



(50) 



Pi P2 1 + PlPl + P2P2 



Here the relation (|4"0"|) has been used. Eq. (|49[) defines a regular Ricmann-Hilbcrt problem (i.e. without zeros). 
Next we take the variation of the Riemann-Hilbert problem (|49| , and get 



SF+ =5F-G + F-5G 1 
Utilizing Eq. (|49|) . we can rewrite the above equation as 



sf.f: 



5F-F1 



f_5Gf; 



(51) 



(52) 



which defines yet another regular Riemann-Hilbert problem for SFF^ 1 . Unlike the previous Riemann-Hilbert problems 
(|33[) and (|49[) which were in matrix product forms, the present Riemann-Hilbert problem (|52[) can be explicitly solved. 
Using the Plcmelj formula, the general solution of this Ricmann-Hilbcrt problem is 



6FF- 1 (C,x)=Mx) + ^ 
Zwi 



n(£aQ 



where 



n(£;x) = F_(£;x) SGfax) F^^x), 



Now we consider the large-£ asymptotics of this solution. Notice from expansions (j3"6")l , (f3"Tj) and relation 

C -> 00, 



(53) 

(54) 
that as 



u(x) 
~2lC~ 



1 + O(i) 
O(I) l + 0(i) ^ 



Li* (x) 

2iQ 



# l + 0(i) 



5F+(C;x) 



O(i) O(I) 



o(i) o(i) ^ 



5u* {x) Su(x) ^/ 1 > 



2/c," 



(55) 



In addition, it is easy to see that 



f-C 



1 

C7-c 



n«;i)df, 



£ — > 00. 



(56) 



When these large-£ expansions are substituted into Eq. (|5"3")) . at 0(1), wc find that A (a;) = 0. At 0(( x ), we get 



Su(x) 



Ui 3 (£]x)d^ Su*(x) 



(57) 
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Now we calculate the elements II13 and II31. From Eq. (|50|) . we see that 



8p x 
5p 2 

Spi 8p 2 PlSpi + p\8pi + p 2 5p2 + P2&P2 



(58) 



(59) 



(60) 



This relation, together with Eq. (|57|) . shows that the variation of the potential Su can be expanded into quadratic 
combinations between Jost solutions <I> and adjoint Jost solutions 3>. This relation is generic in integrable systems. 

Now we calculate explicit expressions for Su. Inserting into (|57[) and recalling our notations (pM|) . we readily 
find that 



Then in view of the F± definitions (|48[) and P± expressions (|T8|) . (|25|) . wc find from (|54[) that 

U(C,x) = J-E(H! + H 2 S)- 1 diag(l,l,s 33 )E- 1 SG E diag(l, 1, s 33 )(H 1 + S~ 1 H 2 )~ 1 E~ 1 JI 1 , 
which simplifies to 

5 pi 
n(C; x)=$\ Sp 2 I 

<S/9 2 



6u 



Here the notations are 





4>kl 




Tpkl 


fe = 


4>k2 


, tpk = 


lpk2 




_ 4>k3 _ 







(031013<5/Ol + 03 1 023 ^2 + 011033^1 + 021033^2) d£. 



4>k=[(t>kX,(t>k2,<j>ks\, 4>k = [i>kl,1pk2,4>k3] , fe=l,2,3. 



Utilizing the symmetry relations (j43|) and (j47|) . the above (Su formula reduces to 

p oo 



(5u 



[(031^13 + 033012) 8 pi + (011033 + 013032) 6 pi] d£. 



(61) 



(62) 



(63) 



This is an important step in our derivation, where symmetry conditions play an important role. Regarding 5u*, its 
formula can be obtained by taking the complex conjugate of the above equation and simplified by using the symmetry 
relations (|46|) . Defining functions 



"31^13 
t>32 013 



"33^12 
^33011 



Zn 



5ll<P33 "t- <Pl3932 
^12033 + 013031 



then the final expression for the variation of the potential (8u,5u*) T is 

[Zi(£;x)6pi(0 + Z 2 (t;;x)6pi(li)]d£. 



5u(x) 


i r°° 


5u*{x) 


T J-oo 



(64) 



(65) 



Notice here that due to the symmetry relations (|4"6"]) , Z 2 {C) is equal to Zj"(C*) with its two components swapped. Also 
notice that Z\ and Z2 are the sum of two terms, where each term is a product of a component of a Jost function and a 
component of an adjoint Jost function. The two-term feature of these functions is caused by the symmetry properties 
of the scattering data and Jost functions, which enable us to combine terms together in the general expansion (|6"Tj) . In 
the massive Thirring model, the counterparts of functions Zip, are also sums of two terms jj]. The two-term structure 
there is not due to symmetry properties, but rather due to the asymptotics of its Jost functions in the spectral plane. 
The feature of each term in Z\ and Z2 being a product between a Jost function and an adjoint Jost function, on the 
other hand, is a generic feature in integrable systems. In previous studies on many integrable equations (such as the 
KdV, NLS, derivative NLS, and massive Thirring equations), it was found that these functions were often "squares" 
or products of Jost functions themselves (thus the name "squared eigenfunctions" ) [l], H, 0, QjJ . This was so simply 
because the spectral operators in those systems were 2x2, for which the adjoint Jost functions (rows of the inverse 
of the Jost-function matrix) are directly proportional to Jost functions themselves. That is not generic however, and 
does not hold for the Sasa-Satsuma equation, or in general for integrable equations whose spectral operator is 3 x 3 
or higher. 
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The derivation of ([65)) for the expansion of (5u,5u*) T is an important result of this subsection. It readily gives the 
variations in the Sasa-Satsuma fields in terms of variations in the initial data. To show this, we first restore the time 
dependence in equation (|65|) . An important fact we need to notice here is that the Jost solutions <f> as defined in (fT3|) 
do not satisfy the time evolution equation of the Lax pair of the Sasa-Satsuma equation. Indeed, this time evolution 
equation of the Lax pair is [23| 



Y t = -4< 3 Ay + V(C, u)Y, 



(66) 



where the matrix V goes to zero as x approaches infinity. Obviously the large- a; asymptotics of $ [see (0 and (|13|) ] 
can not satisfy the above equation as \x\ — ► oo where V vanishes. But this problem can be easily fixed. Defining the 
"time-dependent" Jost functions 



(67) 



then these functions satisfy both parts of the Lax pair, ([2]) and (|66p . The reason they now satisfy the time evolution 
equation (|66p is due to their satisfying the asymptotic time evolution equation (|66p as \x\ — > oo as well as the 
compatibility relation of the Lax pair. Similarly we define the "time-dependent" adjoint Jost functions as 



$(*) = e 4iC 3 At|, 



(68) 



which satisfy both adjoint equations of the Lax pair. Now we use these new Jost functions and adjoint Jost functions 
to replace those in the definitions (|6^|) and get "time-dependent" {Zy^Z-i) functions 



mi 
mi 



mi 
mi 



AD 



mi 
mi 



-"13 ^32 

mi 



(69) 



In view of the relations (|6"T)) and 
ones as 



Another fact we need to notice in conjunction with the time-restored equation of (|65|) is that for the Sasa-Satsuma 
equation, the time evolution of the scattering matrix S is given by [23j 



we see that these "time-dependent" (Zi,Z 2 ) functions are related to the original 
Zf = Z x e^\ zf> = Z 2 e~ 8 « 3t . (70) 



Thus 



S t = -4i?[A,S\. 



»33(*)=0, 4l(*) = 8^ 3 S3lW, 



(71) 



(72) 



where the prime indicates differentiation with respect to t. Then recalling the definition (|42|) of (pi, pi) as well as the 
symmetry property (|10p . we obtain the time evolution of the varied scattering data as 



6 Pl (Z,t)=6p 1 (t,0)e 



(73) 



Now we insert the above time-dependence relations ([70]) and (j73|) into the time- restored expansion (|65|) . and get a 
new expansion 



5u{x, t) 


1 r°° - 


5u*(x, t) 


K ./-oo - 



Zf (£ x, t)5 Pl (f, 0) + Zf {i- x, t)Sp! 0) 



(74) 



The advantages of this expansion are that the Jost functions $^'- ) in [Z{ , Z^ \ satisfy both equations of the Lax pair, 
and the expansion coefficients [Spi,5 P i] are time-independent. 

Since Su and Su* must satisfy the homogeneous linearized Sasa-Satsuma equation, it follows that the functions 
[Z^\ Z^] appearing in the expansion (|74[) must also and precisely satisfy the same linearized Sasa-Satsuma equation. 
A similar fact for several other integrable equations has been noted before by direct verification (see Q for instance). 
Here we will give a simple proof of this fact which holds for any integrable equation where an expansion such as (|74I) 
exists. Suppose the linearized operator of the Sasa-Satsuma equation for [u(x, t), u*(x, t)] T is C. Since the potential 
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u(x,t) satisfies the Sasa-Satsuma equation, and the variation of the potential [6u(x, t), Su* (x, t)] T is infinitcstimal, it 
follows that 



5u(x, t) 
5u*{x, t) 

Inserting Eq. (|74[) into the above equation, we get 

'CZ® (fc x, t)S Pl (£, 0) + Lzf (f ; x, t)5pr (f , 0) 



d£ = 0. 



(75) 



(76) 



Since the initial variations of the scattering data (Spi, Spi)(£, 0) arc arbitrary and linearly independent, it follows that 



CZ®(t;x,t)=CZW(£;x,t)=0 



it). 



(77) 



for any £ <= 



Thus z[^ and Z^' satisfy the linearized Sasa-Satsuma equation and are therefore the squared 



r(t) 



cigenfunctions of this equation. 

Now in the derivation of (|74[) . from (|52p forward, no constraints have been put on the variations of u and u* , except 
for the implied assumption that the scattering data for u and u* and its variations will exist. Whence we expect that 
[5u,Su*] T in (f?4")) can be considered to be arbitrary. If so, then due to the equality which we see in ([73)1. we already 



can expect {Z{ x,£),Z\ s R} to form a complete set in an appropriate functional space (such as L\). 

We remark that the squared cigenfunctions given here are all linearly independent. That can be easily verified upon 
recognizing that from @ , one may readily construct an integro-diffcrential eigenvalue problem, for which these squared 
cigenfunctions arc the true cigenfunctions. Then in the usual manner it follows that these squared cigenfunctions 
are linearly independent. Similar for the adjoint squared eigenfunctions. The linear independence of these squared 
eigenfunctions can also be easily established after we have obtained the dual relations of (j?4"|) which give variations of 
scattering data due to arbitrary variations of the potentials (see next section). 

The squared eigenfunctions [z[*\ Z^] have nice analytic properties. Indeed, recalling the analytic properties of 

Jost functions discussed in Sec. 2, we sec that z[ t} (C)e~ 2 ' Kx is analytic in C_, while Z^\Qe 2iix is analytic in C+. 
These analytic properties will be essential when we extend our results to the general case where S33 and S33 have 
zeros (see end of this section). 

It should be noted that in the expressions (j6"9")l for [z[ , z!f ], if 03 is replaced by other Jost functions (such as 

cf>\ ), or if ^f 1 is replaced by other adjoint Jost functions (such as 1^2 )> the resulting functions would still satisfy 
the linearized Sasa-Satsuma equation and are thus also squared eigenfunctions of this equation. These facts can be 
verified directly by inserting such functions into the linearized Sasa-Satsuma equation and noticing that $W and 
satisfy the Lax pair @, (|6"6"|) and the adjoint Lax pair respectively. Similar results also hold for other integrable 
equations. However, these other squared eigenfunctions in general do not have nice analytic properties in C± nor are 
they linearly independent of the set (f6"9|) . 



B. Variations of the scattering data and adjoint squared eigenfunctions 

In this subsection, we calculate variations of the scattering data which occur due to arbitrary variations in the 
potentials. These formulae, together with Eq. (|6"5|) or (fT4")) . will give the "adjoint" form of the squared cigenfunctions 
as well as their closure relation. 

We start with the spectral equation © for Jost functions <&, or specifically, 

$ x = -<A$ + Q*. (78) 

Taking the variation to this equation, we get 

5^ x = -i(AS$ + Q5$ + SQ<I>. (79) 

Recalling that $ — > e~ z ^ Ax as x — > —00 [see Eqs. ([7]) and (fT5)) ], thus — > as x — > —00. As a result, the solution of 
the inhomogeneous equation (|79p can be found by the method of variation of parameters as 



**(C;a:) = $(C;aO / ^- 1 (C,y)SQ{y)^(C,y)dy. 

J —OO 



(80) 
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Here 



SQ = 



Su 
Su* 
-Su* ~Su 



(81) 



is the variation of the potential. Now we take the limit of x — > oo in the above equation. Recalling $ = if?S and the 
asymptotics of \& — > E as x — > oo, we see that $ — > <5<f> — > £^55 as a; — > oo. Thus in this limit, the above equation 
becomes 



/oo 
<^- 1 ^;x)SQ(x)m^dx, £e 
-oc 

Noticing that S^ -1 = = ^, the above equation can be rewritten as 

poo 

6S(£)= V(£;x)SQ(x)<f>(£;x)dx, £ e R. 



(82) 



(83) 



This formula gives variations of scattering coefficients Ssij via the variation of the potential SQ. From it, the variation 
of the scattering data Spi can be readily found. This Sp\ formula contains both Jost functions $ and adjoint Jost 
functions When we further express $ in terms of iff through the relation (fT5|) (so that the boundary conditions of 
the Jost functions involved are all set at x = +oo), we obtain the expression for Sp\ as 



If 00 - - - - 

Spi = — / [Su (l/>3l¥>3 - + Su* (^32^3 - ^33^l)] dx, 

s 33 J-00 



where the column vector ip is defined as 



(84) 



(85) 



The variation Sp\ can be derived by taking the complex conjugate of the above equation and utilizing the symmetry 
relations (|43|) and (|46|) . Defining functions 



<P = [<Pl, <P2, Va] = $22^1 ~ Sl2^2, 



(86) 



and inner products 



2 



"031^3 - ^33<£>2 
^32^3 - ^33^1 



Oo 



^32^3 - ^33 
^31^3 - "033^2 



(87) 



(f,9)= / f T (x)g(x)da 



then the final expressions for variations of the scattering data Sp\ and Sp\ are 

1 



4.(0 



fii(£;aO, 



Su(x) 
Su*(x) 



(89) 



Spite) 



^33 



(0 



shte;*), 



Su{x) 
5u*{x) 



(90) 



where 

The above two equations are associated with the expansion (|65[) of the previous section, where time is frozen in 
both cases. A small problem with these equations is that the Jost functions and adjoint Jost functions appearing 
in the definitions (|57|) of [Qi , do not satisfy the time evolution equation (|6"6"|) of the Lax pair (see the previous 
section). To fix this problem, we repeat the practice of the previous section and introduce "time-dependent" versions 
of the functions [f2i,f2 2 ] as 



(/) 



4U ] 
4U ] 



4hf 
4U ] 



4U ] 
4U ] 



(91) 
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where (£, x, t) and ip^ (^, x, t) are defined as 

$W and i> (t) have been defined in ([HT |) -([68 ]) . and and are defined as 



(92) 



(93) 



Notice from (f7T|) and the symmetry condition (|10| that (s\2, S22, S21, S22) are time-independent. Similar to what we 
have done in the previous section, we can show that ip^ and in the definition (f9~Tj) are Jost functions satisfying 
both the Lax pair J2]) and ([66]) , while functions yiW and $W satisfy both the adjoint Lax pair. In addition, [fij* , fij ] 
are related to [f2i, ^2] as 



n 



Inserting this relation and (|73|) into the time-restored equations ([89]) and ([90|) . we get the new relations 

1 



5u(x, t) 
5u*{x, t) 



(94) 



(95) 



'33 



(O 



6u(x, t) 
5u*(x, t) 



(96) 



These two new relations are associated with the new expansion (|74[) in the previous section. 

The functions (fi^fi^) appearing in the above variations-of-scattering-data formulae (f95"]l - (j96"|) are precisely ad- 
joint squared eigenfunctions satisfying the adjoint linearized Sasa-Satsuma equation. Analogous facts for several other 
integrable equations have been noted before by direct verification (see @). Below we will present a general proof of 
this fact which will hold for those integrable equations where we have relations similar to l|74[). ([95]) and (|Q6|) . Along 
the way, we will also obtain the closure relation, the orthogonality relations and inner products between the squared 
eigenfunctions and the adjoint squared eigenfunctions. 

First, we insert Eqs. (|9l)|) - (f9"6"|) into the expansion (|T4"|) . Exchanging the order of integration, we get 



Su(x, t) 


2 PQO pQO 


Su*(x, t) 


J — oo J — oo . 



1 



, 3.'i 



(0 



z[ t \^x,t)n { * )T (e 1 x\t) + 



1 



Z®(t;x,t)0!£>*(Z;x',t) 



'33 



(0 



Su(x' , t) 
Su*(x',t) 



dx . 

(97) 

Since we are considering that {5u, Su*) are arbitrary localized functions (in an appropriate functional space such as 
Li), in order for the above equation to hold, we must have 



l —-Z { * ] (£; x, t)n { * )T (e, x', t) + -^Zjpfc x, t)ti 2 t)T (£; x' , t) 



\33 



te) 



'33 



(0' 



d£ = d(x-x')l, 



(98) 



which is the closure relation (discrete-spectrum contributions are absent here due to our assumption of S33 , S33 having 
no zeros). Here S(x) is the dirac delta function. This closure relation has the usually symmetry in that not only 
do the squared eigenfunctions {z[^te; x, i), Z% (£; x, t), £ £ R} form a complete set, but also the adjoint functions 

{fti (£; x, t), CI2 (£; x, t), £ £ M} form a complete set as well. 

To obtain the inner products and orthogonality relations between these functions, we reverse the above and insert 
the expansion ([74]) into Eqs. (|9"5j) - (|9"6"|) . Exchanging the order of integration, we get 



*Pi(f,0) 
*Pi(f,0) 



Spi^, o) 



(99) 



where 



Mte,e,t) = 



.tt) 



nW(C;x I t),zW(r;x I t)) 



nW(C;x I t),zW(^;x I t))\ 



n i * ) te;x,t),zi t) te';x,t)) / 



(100) 
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Since the relation (|99|) is taken to be valid for arbitrary functions of Sp±(£, 0) and <5pi(£, 0), thus 

M{t,?,t) = 18(£-?). (101) 

Consequently we get the inner products and orthogonality relations between squared eigenfunctions [z[*\ zf ] and 
functions [Cl\ , Oj ] as 

(of (£; x, t),Z®g;x,t)) = -7r4,(0<J(£ " O, (102) 
(of & x, t),Z?\t';x,t)) =-nsis(Zm-a (103) 
(Of (£; t), zf (£'; *,*)) = (Of (f ; x, t), Z® (£'; x, i)) = (104) 

for any £, f € K. 

To show that (Oi, O2) are adjoint squared eigenfunctions which satisfy the adjoint linearized equation of the Sasa- 
Satsuma equation, we first separate the temporal derivatives from the spatial ones in the linearized operator C as 

£ = ld t +L, (105) 

where L only involves spatial derivatives. The adjoint linearized operator is then 

C A = -W t +L A , (106) 

where L A is the adjoint operator of L. Now we take the inner product between the equation Cz[ (£';x, t) = and 
function of (£;x,t). When Eq. (|105[) is inserted into it, we find that 

(LZ^Xt 1 ), of (0) - (zPtf), 9 t of (0) + dt(zP (O, fif (0> = o- (107) 

Here and immediately below, the (x,t) dependence of Z± \ an d ^f ^ s suppressed for notational simplicity. Recall 
that S33 is time-independent for the Sasa-Satsuma equation, thus from the inner-product equation (|102[) we see that 
9t(zf (£'), of (£)) = 0. Using integration by parts, the first term in the above equation can be rewritten as 

{Lzf (£'),0f (0) - (Zfta^f (0> + ^',0lx=- i! (108) 

where function ^; a;, t) contains terms which are generated during integration by parts. These terms are quadratic 

combinations of (£'), of (£) and their spatial derivatives. When x = ±L — > ±00, each term has the form 
/(£')<7(£)exp(±i£' L ± i£L ± 4i£' 3 i ± 4i^ 3 t) due to the large-x asymptotics of Jost functions <& and '5. Here /(£') and 
g(£) are related to scattering coefficients and are continuous functions. Then in the sense of generalized functions, 
the last term in Eq. (|108[) is zero due to the Ricmann-Lcbcsgue lemma. Inserting the above results into Eq. (|107p . 
we find that 

(zl t Xe),C A tl®{t)) = 0, ( 6 R, (109) 

i.e. £ A of (£) is orthogonal to all z[ (£'). By taking the inner product between the equation £Z.f (£') = and 
function of (£) and doing similar calculations, we readily find that 

(zf (£'),£ A of (0> = o, ?,SeR. (110) 

In other words, £ A of (£) is also orthogonal to all zf (£')■ Since {^f (£'), #f (£')>£' e R l forms a complete set, 
Eqs. (|109f> - (|l 10[) dictate that £^of has to be zero for any £, x and t, i.e. 

£ A of (C;x,i) = 0, £eR. (ill) 

Thus Of (£) is an adjoint squared eigenfunction that satisfies the adjoint linearized Sasa-Satsuma equation for every 
( £ 1. Similarly, it can be shown that of (£) is also an adjoint squared eigenfunction for every £ 6 R. In short, 
{Of (£))Of £ R} is a complete set of adjoint squared eigenfunctions. 
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Like squared eigenfunctions (z[^ , Z 2 ^ ) , these adjoint squared eigenfunctions (fi^ , f2 2 *^ ) are a ^ so quadratic com- 
binations of Jost functions and adjoint Jost functions. In addition, they have nice analytic properties as well. To 
see this latter fact, notice from the definition ([55)1 that pe^ x is analytic in C_ since 821,8221 V*i e an d ^2 e are 
analytic in C_. Then recalling that Tp^e 1 ^ is also analytic in C_, we see that f2^(£)e 2i ^ x is analytic in C_. Similarly, 
Q^(0e- 2iCx is analytic in C+. 

It is worthy to point out that in the above adjoint squared eigenfunctions (f9Tj) . if tp® is replaced by the Jost 
function ip^ (or ip 2 ), or ^ i s replaced by the adjoint Jost function (or ip 2 ^), the resulting functions would 
still satisfy the adjoint linearized Sasa-Satsuma equation. Thus one may be tempted to simply take ipW = ipjfi, 

(pi*) — (k = 1 or 2) rather than (|9"!2)) as adjoint squared eigenfunctions. The problem with these "simpler" adjoint 
squared eigenfunctions is that their inner products with squared eigenfunctions (|69p are not what they are supposed 
to be [see (|102[) - (|104[) ]. neither does one get the closure relation (|98|) . Thus the adjoint squared eigenfunctions (|9Tj) 
coming from our systematic procedure of variation calculations are the correct ones to use. 



C. Extension to the general case 



In the previous subsections, the squared eigenfunctions and their closure relation were established under the as- 
sumption that S33 and S33 have no zeros, i.e. the spectral equation ([2]) has no discrete eigenvalues. In this subsection, 
we extend those results to the general case where S33 and S33 have zeros in their respective planes of analyticity C± . 
Due to the symmetry properties (|10[) and (fT2")l . we have 

533(0=4^(0, «33(0 = *33(-C*)- (112) 

Thus if Q G C_ is a zero of S33, i.e. s^lQ) = 0, then —Q G C_ is also a zero of S33, and {£*, — (j} G C+ are both 
zeros of S33. This means that zeros of S33 and S33 always appear in quadruples. Suppose all the zeros of S33 and S33 
are Q G C_ and Q G C+, j = 1, 2N, respectively. For simplicity, we also assume that all zeros are simple. 

In this general case with zeros, squared eigenfunctions z[*\ Z^ clearly still satisfy the linearized Sasa-Satsuma 
equation, and adjoint squared eigenfunctions 0,^,0,^ still satisfy the adjoint linearized Sasa-Satsuma equation — 
facts which will not change when S33 and S33 have zeros. The main difference from the previous no-zero case is 
that, the sets of (continuous) squared eigenfunctions {z[ (£;x),Z 2 (£;^)>£ G K} and adjoint squared eigenfunctions 
{Q.^\^;x),Ct 2 {£;x),£ G R} are no longer complete, i.e. the closure relation ([98]) docs not hold any more, and 
contributions from the discrete spectrum must be included now. To account for discrete-spectrum contributions, one 
could proceed by adding variations in the discrete scattering data in the above derivations. But a much easier way 
is to simply pick up the pole contributions to the integrals in the closure relation (|98|) , as we will do below. Recall 
from the previous subsections that z[*\c x, t)e~ 2l ^ x and Q^(C; x, t)e 2% ^ x are analytic in C_, while Z 2 {C,; x, t)e 2l< * x 
and fl 2 t \(;x,t)e~ 2i<: ' x are analytic in C+. In addition, from the large-£ asymptotics (|3"Tj) - (f3"!?)) . ([55]) . and symmetry 
relations (|46[) . we easily find that 

Z 1 t) (C;x,t)^e 2 ^ +8 < 3t (0,l) T , Z i 2 t) (C,x,t)^e- 2 ^-^ 3t (l,0) T , ( 00, (113) 



oW(C;x,t)^-e- 2 ^- 8i ^(0,l) T , fi< t >(C;*,t)^-e 2 * a+8 « 8t (l,0) T , C 00. (114) 

Thus, 



C~ *33 



\-Z[ t) (C,x,t)nf )T {C,x',t)dC = - f e 2i « x ~ x 'Uc diag(0, 1), (115) 



4s(0' 



where the integration path C~ is the lower semi-circle of infinite radius in anti-closewise direction. Since the function 
e 2i({x-x ) j n t ne rigi^ hand gide f the above equation is analytic, its integration path can be brought up to the real 
axis (el. In the sense of generalized functions, that integral is equal to ttS(x — x'), thus 

\prZ[* \C, x, t)n[ t)T (C x', t)d( = -nS{x - x')diag(0, 1). (116) 



C~ *33 



Doing the same calculation for the integral of Z 2 ^ fl 2 ^ T / s 2 3 along the upper semi-circle of infinite radius C + in 
closewise direction and combining it with the above equation, we get 



1 



^zW(C;z,t)fi£ t)T (C;^,tK-- / ^ 7F rZi\ix,t)tl® T (Crf,t)dC = 6{x-x>)l. (117) 

C- «33(C) * JC+ s 33(0 
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Now we evaluate the two integrals in the above equation by bringing down the integration paths to the real axis and 
picking up pole contributions by the residue theorem. Recalling that zeros of S33 and S33 are simple, the poles in 
the above integrals are second-order. After simple calculations, we get the closure relation for the general case (with 
discrete-spectrum contributions) as: 



\»Z®{ti; x, t)rt* )T (Z; x', t) + -±-Z®& x, t)^ )T '(£ x', t) 



2N 



2 1 



^ s 33(0) 



Z®(&;x, t)Q^ T (Q;x', t) + x, t)Q^ )T (Q; x' , t) 

z^^x, t)o$ )T (( r> x', t) + ^§M0; x, m¥ T {Q ]X ', t) 



(/) 



2 A 



T.7. 



2 1 



S(x - x')t, 



(118) 



where 



s 33(0) 



(119) 



Q^iCr^t) = ^(Cr,x,t) - $M^\ Cf ,x,t). (120) 

which have been termed "derivative states" [l], Q due to the differentiation with respect to £. Clearly, Z-f (Cj), 

dz[^ /dC, (Cj), ^2*^(0) i an d dZ^ /d( (Q) satisfy the linearized Sasa-Satsuma equation and are thus discrete squared 

eigenfunctions, while (0 ) > ®i (0)i ^2 an< ^ ®2 (0) satisfy the adjoint linearized Sasa-Satsuma equation 

and are thus discrete adjoint squared eigenfunctions. 

In soliton perturbation theories, explicit expressions for squared eigenfunctions under soliton potentials are needed. 
For such potentials, S31 = S32 = S13 = S23 = 0, thus G = 1 in the Riemann-Hilbcrt problem ([551) . Solutions of this 
Ricmann-Hilbcrt problem have been solved completely for any number of zeros and any orders of their algebraic and 
geometric multiplicities (25| . In the simplest case where all zeros of the solitons are simple, solutions P± can be found 
in Eq. (76) of [25| (setting A 1 = r^' = 1). One only needs to keep in mind here that the zeros of S33 and S33 always 
appear in quadruples due to the symmetries (|112[) . Once solutions P± are obtained, the solitons u(x,t) follow from 
Eq. (|37p . Taking the large-x asymptotics of P± and utilizing the symmetry conditions (|10p . the full scattering matrix 
S can be readily obtained. Then together with P±, explicit expressions for all the Jost functions <f>, ^ as well as 
squared eigenfunctions are then derived. 



IV. SUMMARY AND DISCUSSION 



In this paper, squared eigenfunctions were derived for the Sasa-Satsuma equation. It was shown that these squared 
eigenfunctions arc sums of two terms, where each term is a product of a Jost function and an adjoint Jost function. 
The procedure of this derivation consists of two steps: one is to calculate the variations of the potentials via variations 
of the scattering data by the Riemann-Hilbert method. The other one is to calculate variations of the scattering data 
via variations of the potentials through elementary calculations. It was proved that the functions appearing in these 
variation relations are precisely the squared eigenfunctions and adjoint squared eigenfunctions satisfying respectively 
the linearized equation and the adjoint linearized equation of the Sasa-Satsuma system. More importantly, our 
proof is quite general and generally holds for other integrablc equations. Since the spectral operator of the Sasa- 
Satsuma equation is 3 x 3, while its linearized operator is 2 x 2, wc demonstrated how the two-component squared 
eigenfunctions arc built from only selected components of the three-component Jost functions, and we have seen that 
symmetry properties of Jost functions and the scattering data play an important role here. 

The derivation used in this paper for squared eigenfunctions is a universal technique. Our main contributions to 
this method are triple-fold. First, we showed that for a general integrable system, if one can obtain the variational 
relations between the potentials and the scattering data in terms of the components of the Jost functions and their 
adjoints, then from these variational relations one can immediately identify the squared eigenfunctions and adjoint 
squared eigenfunctions which satisfy respectively the linearized integrablc equation and its adjoint equation. Second, 
we showed that from these variational relations, one can read off the values of nonzero inner products between 
squared eigenfunctions and adjoint squared eigenfunctions. Thus no longer is it necessary to calculate these inner 
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products from the Wronskian relations and the asymptotics of the Jost functions. Thirdly, we showed that from these 
variational relations, one can immediately obtain the closure relation of squared eigenfunctions and adjoint squared 
cigenfunctions. After squared eigenfunctions of an integrable equation are obtained, one then can readily derive the 
recursion operator for that integrable equation, since the squared eigenfunctions are eigenfunctions of the recursion 
operator @, 0, HO, 01 . 

An important remark about this derivation of squared cigenfunctions is that it uses primarily the spectral operator 
of the Lax pair. Indeed, the key variation relations (f6"5")) . and (|5D|) were obtained exclusively from the spectral 
operator ||2J), while the "time-dependent" versions of these relations (|74|) . ([55)) and were obtained by using 
the asymptotic form of the time-evolution equation in the Lax pair as \x\ — > oo. This means that all integrable 
equations with the same spectral operator would share the same squared eigenfunctions (except an exponential-in- 
time factor which is equation-dependent). For instance, a whole hierarchy of integrable equations would possess the 
"same" squared eigenfunctions, since the spectral operators of a hierarchy are the same. This readily reproduces 
earlier results in (Tol . [ill . [l9| , where such results were obtained by the commutability relations between the linearized 
operators and the recursion operator of a hierarchy. 

Since the derivation in this paper uses primarily the spectral operator of the Lax pair, if two integrable equations 
have similar spectral operators, derivations of their squared eigenfunctions would be similar. For instance, the spectral 
operator of the Manakov equations is very similar to that of the Sasa-Satsuma equation (except that the potential 
term Q in the Sasa-Satsuma spectral operator ^ possesses one more symmetry). Making very minor modifications 
to the calculations for the Sasa-Satsuma equation in the text, we can obtain squared cigenfunctions for the Manakov 
equations. This will be demonstrated in the Appendix. 

With the squared cigenfunctions obtained for the Sasa-Satsuma equation, a perturbation theory for Sasa-Satsuma 
solitons can now be developed. This problem lies outside the scope of the present article, and will be left for future 
studies. 
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Appendix: Squared eigenfunctions for the Manakov equation 



The Manakov equations are [26J] 

iu t + u xx + 2{\u\ 2 + \v\ 2 )u = Q, (121) 



ivt + v xx + 2(|u| 2 + \v\ 2 )v = 0. (122) 
The spectral operator for the Manakov equations is Eq. © with 









u 








V 


-u* 


— V* 






Q = v . (123) 
y -u* -v* J 

Below we will use the same notations on Jost functions and the scattering matrix as in the main text. Since Q 
above is also anti-Hermitian, Jost functions and the scattering matrix satisfy the involution properties © and (jTOJ) 
as well. The main difference between the Sasa-Satsuma equation and the Manakov equations is that, the potential 
Q above for the Manakov equations does not possess the additional symmetry ([4J, thus the Jost functions and the 
scattering matrix do not possess the symmetries (| 1 1 1) and (|12p . In addition, the asymptotics ([55]) also needs minor 
modification. It is noted that the spectral operator of the Manakov equations is 3 x 3, while the linearized operator 
of these equations is 4 x 4. So one needs to build four-component squared eigenfunctions from three-component Jost 
functions. This contrasts the Sasa-Satsuma equation where one builds two-component squared eigenfunctions from 
three-component Jost functions. 

Repeating the same calculations as in the main text with the above minor modifications kept in mind, we still get 
Eq. (f5"T)) for the variations (Su, Su*), while the expressions for variations (Sv,Sv*) are Eq. (|5"T|) with II13 and II31 
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replaced by II23 and II32. The expression for matrix II is still given by Eq. ([60]) . Thus variations of the potentials 
via variations of the scattering data are found to be 



[5u, 8v, Su* , Sv 



(124) 



where 
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_ 0131 
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*32 



(125) 



Variations of the scattering matrix SS are still given by Eq. (|83|). with modified in view of the form of (|123|) . 
Then defining vectors 



Pi = {Plj) = S22-01 - Sl2^2 

we find that variations of the scattering data are 



Pi = (Plj) = S 2 2i'l - S21-02, ¥2 = (lf2j) = -Sl2"01 + ^11^2, 

(p~2j) = S2li>l + Sn^2, 



f2 



i 



5 (0 



tti(£;x), 



8u(x) 
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5u*{x) 
Sv*(x) 
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Su*(x) 
Sv*(x) 



(126) 
(127) 

(128) 



where 
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, tt 4 = 
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For Manakov equations, the time evolution operator of the Lax pair is 

Y t = -2iQ 2 kY + V(C, u)Y, 



(129) 



(130) 



(131) 



where V goes to zero as x approaches infinity. Introducing "time-dependent" Jost functions and adjoint Jost functions 
such as 



M 



$(*) = e 2 < 2At $ 



(132) 



and use them to replace the original Jost functions and adjoint Jost functions in the Zk and flk definitions Q125P and 
(|130[) . the resulting functions z[ 2 3 4 and fi^a 3 4 would be squared eigenfunctions and adjoint squared eigenfunctions 
of the Manakov equations. Equivalently, {Zie~ 4ll > i^e -4 '^ 'ji^e 4 ^ *,Z4e 4ll > *} are squared eigenfunctions, and 
{flie 4 ^ *,f22e 4i< ' *, f^e -4 ^"', f^e -4 ^"*} adjoint squared eigenfunctions, of the Manakov equations. 

It is noted that squared eigenfunctions for the Manakov equations have been given in [18[ under different notations 
and derivations. Our expressions above are more explicit and easier to use. 
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